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We describe a search for renormalization group fixed points which control a second-order 
quantum phase transition between a d^2_yi superconductor and some other supercon- 
ducting ground state. Only a few candidate fixed points are found. In the finite tempera- 
ture (T) quantum-critical region of some of these fixed points, the fermion quasiparticle 
lifetime is very short and the spectral function has an energy width of order ksT near 
the Fermi points. Under the same conditions, the thermal conductivity is infinite in the 
scaling limit. We thus provide simple, explicit, examples of quantum theories in two 
dimensions for which a purely fermionic quasiparticle description of transport is badly 
violated. 

1. Introduction 

The quasiparticle excitations of the d^^ _y2-wa?fe high temperature superconductors 
have been subjected to intense scrutiny in the past few years. An especiahy impor- 
tant test of our understanding of the underlying physics is whether quasiparticle 
relaxation processes measured by different experimental probes can be reconciled 
with each other. A striking dichotomy appears to have emerged recently in such 
a context. Photoemission experiments^ indicate that the nodal quasiparticles have 
very short lifetimes in the superconducting state, with their spectral functions hav- 
ing linewidths of order ksT . In contrast, transport experiments, especially mea- 
surements of the thermal conductivity, when interpreted in a simple quasiparticle 
model, indicate far larger quasiparticle lifetimes in the superconductor.^ 

Motivated primarily by the photoemission experiments,^ we recently proposed 
a scenario^'^ in which the large linewidths in the fermion spectral functions are 
explained by fluctuations near a quantum critical point between the dx^~y^ su- 
perconductor and some other superconducting state X (see Fig [ij . In this paper 
we will review and extend earlier arguments which classify various possibilities for 
the state X . We will provide the details of a renormalization group (RG) analysis 
which shows that only a small number of the candidates for X are associated with 
a RG fixed point which describes a second-order phase boundary between X and 
the d^2_y2 superconductor that can be generically crossed by varying a single pa- 
rameter r; we will identify the subset of these fixed points which lead to fermion 
spectral linewidths of order ksT. We will also initiate a discussion of the transport 
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Fig. 1. Finite-temperature phase diagram^ of the d-wave superconductor close to a quantum- 
critical point. Superconductivity is present for T < Tc. The long-range order associated with the 
state X vanishes for T > Tx , but fluctuations of this order provide anomalous damping of the nodal 
quasiparticles in the quantum-critical region. The tuning parameter r is some coupling constant 
in the Hamiltoiiiaii, with the suitable choice depending upon the identity of the state X. It is 
possible, although not necessary, that increasing r corresponds to increasing doping concentration, 
5. The critical point r = rc is not required to be in the experimentally accessible parameter regime, 
just not too far from the d^2_„2 superconductor. 



properties of these fixed points, and argue that they offer attractive possibihties for 
explaining the transport experiments as well. 

Wc begin with a brief review of studies of quantum phase transitions in the 
cuprate superconductors. The subject was initiated in the work of Chakravarty et 
al.,^ who presented a field- theoretic study of the destruction of Neel order in insu- 
lating antifcrromagnets, but focused mainly on thermal fluctuations above a Neel 
ordered state. Subsequent work examined the nature of the paramagnetic ground 
state in the insulator,^ and of its quantum-critical point to the Neel state. ''■^ In par- 
ticular, it was proposed^ that destroying the Neel order by adding a finite density 
of mobile charge carriers also led a quantum-critical point in the same universality 
class as in an insulating antiferromagnet, with dynamic critical exponent z = 1. 
This scenario had strong consequences for magnetic experiments, and for the man- 
ner in which a spin 'pseudo-gap' would develop as T was lowered, and these appear 
to be consistent with observations: crossovers in NMR relaxation rates, ^ uniform 
susceptibility,^ and dynamic neutron scattering^'' all indicate that z = 1. More- 
over, a further consequence'"'^^ of such a scenario was that the paramagnetic state 
should have a stable S = 1 'resonant' spin excitation near the antiferromagnetic 
wavevector, and this is also borne out by numerous neutron scattering studies. 

If we accept that the mobile charge carriers have a superconducting ground 
state (in particular, a rf-wave superconductor), then the arguments for the common 
universality of the magnetic quantum critical point in insulating and doped anti- 
fcrromagnets can be sharpened. For both cases, it is clear that the order parameter 
is a 3-component real field, Na, which measures the amplitude of the local anti- 
ferromagnetic order. In the paramagnetic state, will fluctuate about = 0, 
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and indeed, the triplet 'resonance' modes just mentioned are the 3 normal mode 
oscillations of Na ■ The transition to the state with magnetic long-range order (we 
assume that the charge sector is superconducting on both sides of the transition) is 
described by the condensation of N^, and the theory for the quantum-critical point 
will depend upon whether the Na couple efficiently to other low-energy excitations, 
not directly associated with the magnetic transition. In a d-wave superconductor, 
the important candidates for these low energy excitations are the fermionic 5 = 1/2 
Bogoliubov quasiparticles (one can also consider fluctuations in the phase of the su- 
perconducting order parameter, but we will argue below that such a coupling can 
be safely neglected). Momentum conservation now plays a key role: fluctuations of 
Na occur primarily at a finite wavevector Q (in the present situation, this is the 
wavevector at which the magnetic order appears), and the fermions will be scat- 
tered by this the wavevector. If Q does not equal the separation between two nodal 
points of the d-wave superconductor [the nodes are the locations in the Brillouin 
zone of gapless fermionic excitations, and they are at (iif, zLK) with K = O.SQtt at 
optimal doping], then it is not difficult to show that the fermion scattering serves 
mainly to renormalize the parameters in the effective low energy action for the Na, 
and does not lead to any disruptive low energy damping. '^'^ In such a situation, 
there is no fundamental difference between the magnetic fluctuations in a super- 
conductor and an insulating paramagnet, and both cases have the same theory for 
the quantum critical point to the onset of long-range magnetic order. Conversely, 
if Q = (2K,2K), the coupling to the fermionic quasiparticles is important, and a 
new theory obtains: such a theory was discussed by Balents et al}'^ 

In this paper we will focus on the fermionic quasiparticles rather than the mag- 
netic excitations. We are interested in damping mechanisms for the fermions and 
associated possibilities for the state X in Fig [T] The above discussion on magnetic 
properties suggests a natural possibility for X: a state with co-existing magnetic 
and superconducting order. However, it should also be clear from the discussion 
above that strong damping of the fermionic quasiparticles requires Q — {2K, 2K). 
For the current experimental values, this condition is far from being satisfied, and 
so the magnetic ordering transition is just as in an insulator. This also means that 
the magnetic quantum critical point is not currently a favorable candidate for the 
quantum critical point in Fig [TJ We are therefore led to a search for other possi- 
bilities, and this paper will describe the results of such a search. The new cases we 
will consider do not have magnetic order parameters, and so we are envisaging two 
distinct quantum critical points near the d-wave superconductor: one involving the 
magnetic order parameter which is already known to occur with decreasing doping 
(and which, as discussed above, is in accord with numerous magnetic experiments), 
and another one associated with the state X which may or may not be present along 
the experimentally accessible parameter regime. For completeness, we will also dis- 
cuss the magnetic case with Q = {2K, 2K) in Section 12.21 this is the only case 
which envisages a single quantum phase transition to explain both the magnetic 
experiments and the fermion damping. 
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In Section [21 we present the details of a renormalization group (RG) analysis 
which searches for candidate fixed points describing the quantum phase transition 
between X and the d-wave superconductor: among the many a priori possibilities, 
only a few stable fixed points are found. The structure of the T > single fermion 
Green's functions near these fixed points is discussed in Section O these results can 
be compared with photoemission experiments. Spin, thermal, and charge transport 
properties of all the fixed points are subsequently considered in Section [U a signifi- 
cant result is that the thermal conductivity of the quantum field theories describing 
the vicinities of these fixed points is infinite in the scaling limit. 

2. Renormalization group analysis 

Let us assume that the order parameter associated with X carries total momentum 
Q. We have argued^'^ above that its coupling to the fermions can be relevant only 
if a nesting condition is satisfied: order parameter fiuctuations will scatter fermions 
by a momentum Q, and such scattering events are important only if they occur 
between low energy fermions, i.e., the wavevector Q (or an integer multiple of Q) 
must equal the separation between two nodal points. If such a condition is not 
satisfied, then, as noted above, fermion scattering events can be treated as virtual 
processes which modify the coupling constants in the effective action, but do not 
lead to a fundamental change in the form of the low energy theory. 

Three natural possibilities can satisfy the nesting condition:^^ Q =0, Q = 
{2K,2K), and Q = (2i^, 0), (0, 2i4r). We will consider these possibilities in turn 
in the following subsections. Of these, the first condition can be naturally satisfied 
for a range of parameter values, while the last two require fine-tuning unless, for 
some reason, there is a mode-locking between the values of Q and K. 

2.1. Order parameters with Q = 

We will assume that the Q = order parameter is a spin-singlet fermion bilinear 
(spin triplet condensation at Q = would imply ferromagnetic correlations which 
are unlikely to be present, while order parameters involving higher-order fermion 
correlations are not expected to have a relevant coupling to the fermions^''^). Simple 
group theoretic arguments'* permit a complete classification of such order parame- 
ters. The order parameter for X must be built out of the following correlators (cqa 
annihilates an electron with momentum q and spin a —] , |) 

(''qo'^qa) = 

(Cq|C_q|> = [Ao(cOS(7a; - COaqy) + Bq] 6**^, (1) 

where Aq is the background dx-2_yi pairing which is assumed to be non-zero on 
both sides of the transition, ip is the overall phase of the superconducting order, 
and ^q and -Bq contain the possible order parameters for the state X corresponding 
to condensation in the particle-hole (or excitonic) channel or additional particle- 
particle pairing respectively. It is clear that (p has the usual charge 2 transformation 
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(A), (B) 

fully gapped spectrum (G) 




Fig. 2. Evolution of the nodal points when passing from r > Tc to r < rc at T = for 6 of the 
Q = states X discussed in Section 12.11 There is no change in the position of the nodal points 
for case C, only a change in velocity. 

under the electromagnetic gauge transformation, and so gradients of (p measure flow 
of physical electrical current; because of the non-zero superfluid density associated 
with Aq, (p fluctuations remain non-critical and we will show that they can be 
neglected. It follows that the order parameter Bq, which is in general a complex 
number, carries no electrical charge. Similarly, Aq is also neutral but must be real. 

To classify the order parameters for X, we expand Aq and Bq in terms of 
the basis functions of the irreducible representation of the tetragonal point group 
Ciy. This group has 4 one-dimcnsional representations, which we label as (basis 
functions in parentheses) s (1), dx2_y2 {cosq^ — cosgj,), dxy (sing^; sing^), and g 
(sin sin qy (cos — cos % ) ) , and one 2-dimensional representation p (sin q^ , sin qy ) . 
An analysis of all excitonic order and additional pairings (both real and imaginary) 
in these functions has been carried out,^ and it was found that 7 inequivalent order 
parameters are allowed for X. Of these, the first 6 (A-F) involve a one-dimensional 
representation of C^y , and so the order parameter is Ising-like and represented by a 
real field 0, while the 7th (G) involves the 2-dimcnsional representation and 2 real 
fields (pxj't'y The 7 possibilities for X are 

(A) is pairing: Aq = 0, Bq = i(f) 

(B) idxy pairing: = 0, Bq = i(j)smqx sin qy 

(C) ig pairing: Aq = 0, Bq = i(j>smqx sin (cos q^; — cosgj,) 

(D) s pairing: Aq — 0, Bq — (I) 

(E) dxy excitons: Aq = cf) sin q^ sin qy , Bq = 

(F) dxy pairing: Aq = 0, Bq = cf) sin q^ sin qy 

(G) p excitons: ^q = 4>xisinqx + sinqy) + (j)y{sinqx - sinqy), Bq = (2) 

Apart from case C, these possible orderings for X change the nature of the nodal 
excitations, and these are sketched in Fig [2l Cases A, B open up a gap in the 
fermion spectrum over the entire Brillouin zone (suggesting an energetic preference 
for these cases), case C leaves the positions of the nodal points unchanged but 
changes a velocity in the dispersion relation, while cases D-G break C^y symmetries 
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by moving the nodal points as shown. 

Knowledge of the order parameters in ^ and simple symmetry considerations 
allow us to write down the quantum field theories for the transition between the 
dj.2_y2 superconductor and X. 

The first term in the action, S'lp, is simply that for the low energy fermionic 
excitations in the (1^2 _y2 superconductor. We denote the components of Cqa in the 
vicinity of the four nodal points {±K,±K) by (anti-clockwise) fia, /2a, fsa, fia, 
and introduce the 4-component Nambu spinors ^fi^ = {fia, Sabfli,) and ^'20 = 
if 2a, Sabflf,) wherc Eab = —£ba and £|x — 1 [we will follow the convention of writing 
out spin indices (a, 6) explicitly, while indices in Nambu space will be implicit]. 
Expanding to linear order in gradients from the nodal points, we obtain 



/^2^ ^ 



(3) 



Here ujn is a Matsubara frequency, r" are Pauli matrices which act in the fermionic 
particle-hole space, k^^y measure the wavevector from the nodal points and have 
been rotated by 45 degrees from q^^y co-ordinates, and vp^ va are velocities. 

The second term, describes the effective action for the order parameter alone, 
generated by integrating out high energy fermionic degrees of freedom. For cases 
A-F this is the familiar field theory of an Ising model in 2+1 dimensions 



Sd, — 



Sxdr 



rl 



{OA? 



Uo 

24' 



(4) 



here r is imaginary time, c is a velocity, r tunes the system across the quantum 
critical point, and uq is a quartic self-interaction. For case G, the generalization of 
Sd, is 



d^xdr 



{{drCb^f + (dr(by)^ + cl{d^4>xf + Cl{dy4>^f + cl{d,4>yf (5) 



+ Cl{dy4>yf + e{d,<j,.){dy<j,y) + T {<g + •/'^ ) } + ^ {uo(0^ + ) + 2i;O0^0^} 



The final term in the action, S<^ip couples the bosonic and fermionic degrees of 
freedom. From ([2]) we deduce for A-F that 



d^xdr 



\O0 (^-LAfl^-la + 1'L*^2«'2a) 



(6) 



where Aq is the required linear coupling constant between the order parameter and 
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a fermion bilinear. The matrices Mi, M2 are given by 



(A) A'h 


= tV, M2 = ry 


(B) Ml 


= ry, M2 = -ry 


(C) Ao = 


= 


(D) Ml 


= r^, M2 = 


(E) Ml 


= T^ M2 = -T^ 


(F) Ml 


= r^, M2 = -r^ 



(7) 

Note that there is no non-derivative coupUng between (p and for case C: the 
fermions are essentiahy spectators of the transition for this case, which will not be 
considered further. Finally, for case G, 5**0 generalizes to 



J <fxdT [Ao *i + 0^*^*2) ] . (8) 



We can also consider the coupling between the order parameter cj) and the phase 
of the superconducting order, 93, in ([1]). By symmetry, the simplest allowed coupling 
is (jPiyifY. It is easy to show that such a coupling is irrelevant. So the power- law 
correlations generated by the superflow are not important. 

We now describe our RG analysis of the 6 distinct field theories represented by 
S^i + + S^^ and 5** -I- S'^ -I- Sq,^. The familiar momentum-shell method, in which 
degrees of freedom with momenta between A and A — ciA are successively integrated 
out, fails: the combination of momentum dependent renormalizations at one loop, 
the direction-dependent velocities {vp, vai c . . . ), and the hard momentum cut-off 
generate unphysical non-analytic terms in the effective action. So we will construct 
RG equations by^^ using a soft cut-off at scale A, and by taking a K{d/dK) derivative 
of the renormalized vertices and self energies. We write the bare propagators as 

ipi + VFP2T'' + VAPsr'' 2 / A 2n /on 

and 

= ^Kip^A^) (10) 

where K{y) is some decaying cuf-off function with K{0) = 1; e.g., K{y) = e~y 
is a convenient choice. [The momentum-shell method would correspond to a step- 
function cut-off K(y) = 6(y — 1).] To handle possible anisotropics we use a hybrid 
approach; the space-time integrals over p = {ujn,Px,Py) are written down in D = 
2-1-1 dimensions, they can be split via p = pn (where n is a unit vector) into 
an angular integral / dfln containing all direction dependent information and an 
integral over p — |p| which is then evaluated in D = A — e dimensions. 

We demonstrate this method by calculating the linear-order fermionic self-energy 
for the cases A,B,D-F. The diagram shown in Fig[3}3 evaluates to 
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Fig. 3. Diagrams entering the one-loop RG equations for cases AjBjD^F- FuU lines (dashed) lines 
denote fermionic (bosonic) propagators, filled (open) circles represent the Aq (uq) interaction. The 
constants Ai, ... are obtained by evaluating the diagrams, expanding to lowest non-trivial order 
in the external momentum and taking the A{d/dA) derivative as described in the text, a) and b) 
are the velocity renormalizations whereas c) contains the renormalized couplings. (Diagrams for 
^2 are similar.) 



For S^fj each vertex contains the couphng matrix Mi, the signs S2, S3 — ±1 are 
therefore given by Mir^Afi = S2T^, Mir^Afi — s^t^ . Expanding the above expres- 
sion to hnear order in q gives: 



dp 
8^ 



nl + VpTil + v\n\ 

We have inserted a lower hmit /i in the p integral representing an external momen- 
tum to regularized the infrared divergence. Note that the only quantity entering the 
RG equations is the K{d/dK) derivative of the self-energy. This removes the infrared 
divergence, i.e., we can take the limit A//i 00 in the final expression. Writing the 
cut-off integral involving K in general dimension D we obtain 

d _ 2XlKi f inlqi + S2VFn\q2T^ + ssVAulq^T 



dA *^ ^ (27r)^A4-c J n^^+yj^nl+vlnl ' ^ ' 

/>OC 

K,= dy [-K{y)K'{y) + y[K' {y)\' + yK{y)K"{y)] . 

Jo 

A simple integration by parts for D — A shows that Ki equals unity, independent 
olK{y). 

At this point it is useful to introduce dimensionless coupling constants A and u 
by Ao = XA'/^/Sy,'^ and uq ^uA'/Sd with Sd = f^D/(27r)^ = 2/[r{D/2){ATT)^^^] 
and Qd being the area of a Z3-dimensional unit sphere. From Eq (|12|) we can easily 
read off the values of the Bi : 



B.. = -2K, I ^ ,2^^:;- ^,2,2 (13) 



with si = —1. The other diagrams in Fig [3] are evaluated similarly. 
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The derivation of the RG equations is straightforward and essentially identical 
to the momentum shell method. We construct flow equations for all velocities and 
couplings to one-loop order in the non-linearities A, u. It is convenient to choose 
the time scale such that c = 1 at each step (this introduces a non-trivial dynamical 
critical exponent z). The procedure results in the following RG equations for cases 
A,B,D-F: 

f3{vF) = -vfX' (^B2 - Si + ^^^y^ 
PivA) = -vaX^ (^Bs -B, + 

I3{u) = -eu + ^(v4i + 3A2) - y^D, - 2A\^D2 (14) 

with the renormalization constants as defined in Fig [3] and Ai = 2Ai, A2 = A2 + 
A^. At a fixed point the dynamical critical exponent z and the anomalous field 
dimensions are given by 

z = l + A*2(Ai-A2)/2, 

vb - y^A2 , 

= A*2 [B, - {A, - A2)/2] . (15) 

From Eq (|14p we see that /3{vp) = /3{va) — requires B2 — B3 a.t a non-trivial 
fixed point. The signs of the Bi (fTS]) are determined by Si, and this preempts a fixed 
point for the cases D-F: the structure of Mi^2 leads to S2 = —S3, so that B2 and 
B3 have always different signs. 

The cases A and B have S2 = S3 = —1, and it is easy to show that the fixed 
point equations for vp and va are simultaneously satisfied only for Vp = = 1. 
This implies that the resulting fixed point is Lorentz invariant;^ for _D = 4 the 
renormalization constants are given by Ai^2 — 4:, Bi — 1/2, Ci = —1, Di = —36, 
D2 ~ 2. The RG equations have the infrared stable fixed point'^ with z = 1 and 

Let us briefly discuss the remaining case G. By choosing the time scale we keep 
the velocity ci = 1 fixed. The flow equations for vf-, va, and A have the same form 
as above, but with A12 replaced by 2, since each fermion field couples to only 
one species of bosons. The other RG equations read: 

P{e)^eX^A2, (3(4) ^ clX^As - A2) , 

P^u) = -eu + ^(Ai + 3A2) - ^^-^D^ - 24A4i?2 , 

A^?j 111) iP" 

f3{v) = -ev + — (^1 + 3A2) - - -E2 (17) 
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where Ei and E2 renormalizing v arise from diagrams similar to Fig[5}32. Analysis 
of these equations shows that P{vp) — /3(wa) — Pic^) = is only fulfilled for C2 = 1 
and vp = WA- Furthermore we must have e = at a fixed point, stability requires 
A2 > 0. Explicit evaluation in = 4 shows that the zeros of (3{vf) have A2 < Q 
which proves the non-existence of a fixed point for case G. 

We summarize this subsection by restating the main conclusions. Among the 
Q = order parameters considered here, only cases A and B possess stable RG 
fixed points which can describe a second-order quantum phase transition to the 
state X, and with strong damping of quasiparticle spectral functions. The fixed 
points for both cases are Lorentz invariant, and for future convenience, we display 
this explicitly. We perform the unitary transformation 

(this ensures that the structure of the r matrices in ([3]) is the same for and 
^fj), define ja — — iT^^'j and introduce the Lorentz index ^ — T,x,y. Then 
+ Scf, + S'lJ^ can be written as 



- ^ ^ ' ^2a (18) 



-zAo0(*la4'laT*2a*2a)): (19) 



-0^ + 

2^ 2r 



where 7^ = {—t'^, t^,t^), we have set vp — = c = 1 and the — (+) sign in the 
last term is for case A (B). 



2.2. Order parameters with Q = {2K, 2K) 

One plausible candidate for ordering at such a wavevector is antiferromagnetic order 
at Q = (tt, tt), K — 7r/2; this case was discussed in Section [1] and has been ana- 
lyzed by Balents et aZ.^^ (The generalization to magnetic order at incommensurate 
wavevectors is not diflticult, but we will not consider it because incommensuration 
along the diagonal direction in the Brillouin zone has not been observed in the su- 
perconducting cuprates.) We represent the strength of the antiferromagnetic order 
by a real, three-component field Na (a = x,y,z] for the incommensurate case Na 
is complex). The transition from a d-wave superconductor with (No) = to a state 
X which is a d-wave superconductor with (Na) 7^ is described by the following 
continuum action near the critical point 

S2= J d^x[mad^.r^la + *^2a3p7''*2a + ^(^^^^o)' + + ^{N^J^ 

+ l\oN^eac<Jcb [^laC^-lb - ■^2aC-^2b + H-C] ) , (20) 

where a" are Pauli matrices in spin space, and C — ir^ is a matrix in Nambu space. 
Like Si, S2 has the property of Lorentz invariance {'i'C'i' and are Lorentz 
scalars). Also like Si, the couplings A and u approach^^ fixed point values under 
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the RG, and so the fermion spectral function in the quantum-critical region will 
have an energy width of order ksT. 

For completeness, we mention another order parameter at Q = (tt, tt) which has 
been the focus of some recent discussion: the staggered flux order. The coupling of 
this order to the nodal fermions involves a spatial derivative, and can be shown to 
be irrelevant:'^ so the fermionic quasiparticles are merely spectators at a transition 
involving the onset of staggered flux order and such a theory is not of interest here. 

2.3. Order parameter with Q = (2K, 0), (0, 2K) 

An attractive possibility of an order parameter with such a Q is the charge density 
itself (or 'charge stripe' order). We consider the incommensurate case, in which the 
charge density has the modulation 

5p = Re [$^e'2K^ $^e«2/fyj ^21) 

where ^x,y are complex fields which constitute the charge-stripe order parameter. 
Related order parameters have been considered earlier by Castellani et al,^^ but they 
discussed the onset of charge density wave order in a Fermi liquid (for this case K 
has no meaning, and one considers instead charge density wave order at a generic 
Q which can connect two points on the Fermi surface), whereas we will consider the 
onset in a d-wave superconductor. The theories for the two cases are very different: 
the former has an overdamped spectrum for the order parameter fluctuations^^ and 
(because the theory is not below its upper-critical dimension) does not satisfy the 
strong scaling properties assumed for the fermion spectral functions to be discussed 
in Section [3l which do apply to the latter. We believe it is important to discuss the 
onset of charge order in the true ground state of the doped antiferromagnet, the 
d-wave superconductor, and hence our focus on this case. 

The field theory for a transition from a d-wave superconductor to a state X which 
is a d-wave superconductor co-existing with the modulation pT|) can be deduced'^' 
as in the cases above, and near the critical fixed point it takes the form 

+ Y(l*=.|^ + l*yO +^'o|*-n*!;l' + Ao[$.*2a*la-*yeab*2aC*lb + H.C.] ). (22) 

The RG properties"^ of S3, are very similar to the cases 6*1.2 considered previously: 
u, V, and A approach fixed point values, and there is a universal coupling between 
fermionic and bosonic fiuctuations in the critical region. 

3. Fermion spectral functions 

The remainder of this paper will highlight some important observable properties of 
the field theories 5'i^2,3 in P9l20l22p . We will restrict ourselves to the scaling limit, 
both at T = and T > 0; for some quantities we will see that it is necessary to 
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consider corrections to scaling to obtain a complete picture — we will not enter into 
such issues here and leave them for future work. 

In this section we consider single fermion Green's functions which can be mea- 
sured in photoemission experiments. In the quantum-critical region (Fig [T]) of 5*1,2,3 
these will obey 

where we have set h, ks and all velocities to unity. The scale factor Af is non- 
universal, while the exponent rjf and the complex- valued function are universal 
and depend only upon whether the relevant theory is Si , 5*2 or S3 . Scaling functions 
like $ / have been studied by a variety of approximate methods,^" and we will quote 
some useful limiting forms here (in contrast, in l-t-l dimensions exact results for 
analogous scaling functions are known^"'^^ for all values of their arguments). 
For w, fc > T, (HSl) reduces to^-" 

Gf{k,co) = ACf ^ + k,r^ + Y\^ , (24) 
^ ^[P-(u; + zO)2]i-''^/2' ' 

where C/ is a universal number. Note that the imaginary part of this is non-zero only 
for u) > k, and it decays as for large uj; a similar kinematic constraint and 

large u) tail was noted recently in one-dimensional stripe models for fermion spectral 
functions. The expression (|24|) has a singularity precisely at w = fc: this singularity 
is rounded at any non-zero T > and so ((24|) is not an accurate representation of 
([23|) right at the threshold. Unfortunately, the nature of this rounding is not easy 
to compute in general; it can be computed in an expansion in r/j, and these results 
are in Appendix A of our recent work.'^ 

In the opposite limit, w, fc ^ T, (|23p has a very different form. Now the k and lo 
dependencies are smooth, and all anomalous powers involve only factors of T. The 
following approximate form was suggested"^ in this limit 

where the overall scale of the right hand side defines the value of A, and F f /T is 
a universal number. Estimates of this universal number have been given for S1.3. 
Note that along the diagonal direction in the Brillouin zone {k^ = 0), the result 
psp predicts a simple Lorentzian for the spectral function with an energy width of 
order F^. Of course, the Lorentzian form breaks down for large w, when the spectral 
function crosses over to the slowly decaying tail predicted by (f24|) . 



4. Transport properties 

Unlike the single-particle properties of Section [3l transport properties can require 
determination of composite correlators of both the fermionic and bosonic degrees 
of freedom. In particular, we can have processes in which there is rapid scattering 
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between fermions and bosons, and thus a broad fermion spectral function, but es- 
sentially no degradation of the transport current. In such cases, it is clearly not 
appropriate to interpret transport properties in terms of a single fermionic quasi- 
particle lifetime, and attempts to do so will suggest very long quasiparticlc lifetimes 
very different from those observed in photoemission experiments. 

We describe spin, thermal, and charge transport in the following subsections: 



4.1. Spin transport 

The global SU{2) spin-rotation symmetry of 5*1.2,3 implies that total spin is 
conserved. Consequently, there is an associated Lorentz 3-current, J*" obeying 
9^J^" = 0. For S'1,3 this 3-current has only a fermionic contribution 

= I [*la<fc7M*lb + *2a<b7M*26] (26) 

(this expression agrees with earlier work^^), while for S2 there is also bosonic con- 
tribution^" 

1 
2 

The conservation of J^" implies that the total spin is a constant of motion 



dr 

however, the physical spin current is not conserved: 

d 



/ d^xJ'^"" = 0; (28) 



J d'xJ^% ^ 0. (29) 

The spin conductivity, CTs, is expressed by a Kubo formula involving J^^ (along 
with a 'diamagnetic' contact term for the bosonic fields in S2) and ((29)) implies that 
a ballistic contribution to CTs is not expected. The scaling properties of as can be 
deduced along the lines discussed^° for a purely bosonic theory: the conservation of 
the 3-current implies that the scaling dimension of CTs is protected to take the value 
d — 2 = 0, and so in the quantum critical region we have*'^*^ 

= $s (^) , (30) 

where $s is a universal function of order unity (we have absorbed factors of g/Js 
in the definition of as). In principle, it should not be difficult to extend earlier 
results^" for purely bosonic theories to determine the collisionless {ui ^ T) and 
collision-dominated (w <C T) regimes of as ■ 



4.2. Thermal transport 

This section is based on ideas of T. Senthil (unpublished) in a different context. 

The Lorentz invariance of <5'i,2,3 implies the existence of a energy-momentum 3- 
tensor which is conserved, df^Qf^^u = 0, and which can be chosen to symmetric in 
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the jji, V indices. The energy density is 6t-t, while the energy current is Jf = 9ii-, 
where i = x,y. The conservation of the 3-tensor impUes three constants of 
motion, the total energy and the total momentum: 



Now using the symmetry of and in particular — Qri, we conclude that 
the spatial integral of the energy current Jf is also conserved [contrast this with 
(j29p , in which the spatial integral over the spin current was not conserved] . As the 
thermal conductivity is given by a Kubo formula involving the energy current, we 
conclude that heat transport is ballistic and the thermal conductivity is infinite in 
the scaling limits defined by 5i_2,3- 

4.3. Charge transport 

The above theories describe excitations in superconductors, and a complete de- 
scription of charge transport, and an accounting of the conservation of total charge, 
requires inclusion of the Cooper pairs [alternatively stated, we have to include the 
fluctuations of Lp in ([1])]. However, in a frequency regime where the response is 
dominated by quasiparticle excitations, we can consider the predictions of >S'i^2,3 
alone. The quasiparticle contribution to the electrical current is^^ the 2-vector 
Jf = («'l„^'ia,*L*2a) i(Wi7^«'i,-'^2a7r*2a)- 1^ general there is nothing spe- 
cial about this operator: it will acquire anomalous dimensions at the critical point, 
and this will lead to a contribution to the charge conductivity which has prefactors 
of a non-universal constant and an anomalous power of T . However, for the case 
of Si only, an exceptional circumstance occurs: the action has U{1) symmetries 
corresponding to global changes in the phases oi'^i 2: corresponding conserved 
3-currents 9^\E'ia7M^ia — f^M^2Q7M*2a = Note that these 3-currents do not cor- 
respond to the conservation of physical electrical charge: rather, the t component 
of these 3-currents happen to be equal to the x, y components of the electrical 
current Jf. The conservations of these currents implies that the spatial integral of 
Jf is a constant of motion, and hence the quasiparticle contribution to the charge 
conductivity is a zero- frequency delta function for 5*1. 

5. Conclusion 

This paper has identified quantum field theories at which the nodal quasiparticles of 
a d-wave superconductor undergo strong inelastic scattering. All these theories are 
associated with a quantum critical point between a c?-wave superconductor and some 
other superconducting state X (Fig[T]). We performed an exhaustive search among 
states X associated with a spin-singlet, zero momentum, fermion bilinear order 
parameter, and found that only two candidates possessed a non-trivial quantum 
critical point: X a d^'^-y^ + is or a dj.2_j^2 + id^y superconductor. The quantum 
field theory for these cases is 5*1 in ()19ll . Among order parameters with a non-zero 




(31) 
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momentum, we had to permit a fine-tuning condition that the ordering momentum 
was equal to the spacing between two Fermi points of the d-wave superconductor. 
Two plausible cases of this kind had X possessing co-existing spin-density wave 
order and d-wave superonductivity [described by S2 in ((20)) ] or co-existing charge- 
density wave order and d-wave superconductivity [described by 5*3 in (|22p ] . 

We reviewed the single particle (Section [3]) and transport (Section |4]) properties 
of *5'i^2,3- An important observation was that there is no simple relationship between 
the relaxation times associated with these observables. In particular, strong scatter- 
ing between the fermionic and bosonic modes can lead to very short single particle 
lifetimes and broad spectral functions (Section[3]). At the same time, these scattering 
processes can do little to degrade the transport current: the thermal conductivity 
of 'S'i_2,3 was found to be infinite in the continuum scaling limit. 

Comparison of the above results with experiments requires some understanding 
of the effects of corrections to scaling, and of impurity scattering. Initial steps in 
this direction have been taken, and studies for transport properties are in progress. 

Keeping the above cautions in mind, we offer a possible interpretation of recent 
thermal conductivity measurements.^ A strong enhancement of the thermal Hall 
conductivity, Kxy, is observed as T is lowered, but K^y saturates and decreases 
below T = 28 K. The enhancement suggests large fiuctuations of chiral order, as 
is the case when X is a (1^2 _y2 + idxy superconductor (for T < Tx this state has 
a non-zero K^y even in zero field^^). So we suggest that the experimental system is 
in the vicinity of such a state X but has r > Vc, and the crossover at the dashed 
line in Fig [T] occurs at T « 28 K. A testable implication is that the photoemission 
linewidths of the nodal fermions should be ^ ksT down to 28 K, and then cross 
over to behavior expected in an ordinary BCS o?j,2_j,2 superconductor at lower T. 
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Erratum: Quantum phase transitions in d-wave superconductors 
Phys. Rev. Lett. 85, 4940 (2000) 

Matthias Vojta 

Institut fiir Theoretische Physik, Universitdt zu Koln, Ziilpicher Str, 77, 50937 Koln, Germany 

Ying Zhang 

Goldman Sachs, 1 New York Plaza, New York, NY 10004 
Subir Sachdev 

Department of Physics, Harvard University, Cambridge, MA 02138 

We correct an error in our paper Phys. Rev. 

Lett. 85, 4940 (2000) arXiv:cond-mat/0007170 . Our characterization of the physical 
properties of the superconducting state G was incorrect: it breaks time-reversal symme- 
try, carries spontaneous currents, and possesses Fermi surface pockets. 



In the discussion of case G in our paper, ^ above Eq. (4), the single sentence 
"The state X retains T and the gapless nodal points, but has C^y broken to Z2" is 
incorrect. The state X = G breaks T (time- reversal), and has spontaneous electrical 
currents. For <f>x ^ Q and (f)y = (or vice versa) the currents have the same symmetry 
as those in the state Qu discussed by Simon and Varma.^ Also, as pointed out by 
Berg et al.,^ the nodal quasiparticles do not survive in the superconducting state G, 
but turn into Fermi pockets. The latter conclusion can be verified from the fermion 
spectrum obtained by diagonalizing Eqs. (l) + (5) for constant (l)x,y 

All other sentences and the conclusions in the paper^ remain unchanged. 

Also, in the companion paper, ^ the only error is in the sketch of the fermion 
excitations in Fig. 2 for case G. 

We thank Erez Berg and Cenke Xu for helpful discussions. This research was 
supported by NSF grant DMR-0537077 and DFG SFB 608. 
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